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Abstract 

Two-dimensional disordered quantum antiferromagnets are stud- 
ied by means of a continuum description in which disorder is in- 
troduced by a random distribution of couphngs (spin stiffnesses) in 
the ordered phase of the Nonhnear Sigma Model. Quenched soliton 
(skyrmion) correlation functions are evaluated and used, along with 
the quenched magnetization, to characterize the phase structure of 
the system. When magnetic dilution is exponentially suppressed, the 
introduction of disorder only modifies the subleading terms in the 
large distance behavior of the soliton correlation functions, yielding 
the same skyrmion energy as in the pure case. The system is in a 
"hard" disordered Neel phase similar to the ordered antiferromagntic 
phase occurring in the pure case. Conversely, when magnetic dilu- 
tion is not exponentially suppressed, the large distance behavior of 
the correlation functions is drastically changed. The system exists 
in a new phase in which the energy of quantum skyrmions is equal 
to zero in spite of the existence of a nonvanishing antiferromagnetic 
order parameter. This "soft" disordered Neel phase is characterized 
by universality classes which are determined by the behavior of the 
distribution of random couplings in the small coupling region. The 
possible relation of this phase with spin glasses is briefly discussed. 
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1) Introduction 



The continuum description of two-dimensional quantum antiferromagnets has been 
the object of intense investigation for a long period of time until quite recently 
|jl|, 1^, 1^, 1^. The interest in this kind of description has been mostly enhanced be- 
cause of its successful applications in the case of layered antiferromagnets such as 
the high temperature superconducting cuprates. For these compounds, the undoped 
parent materials can be very well described by the two-dimensional antiferromag- 
netic Heisenberg model on a square lattice [^]. In the continuum limit, this can be 
mapped into the ordered phase of the Nonlinear Sigma Model (NLSM) whose 
single coupling constant ps, the spin stiffness, is directly related to the Heisenberg 
antiferromagnetic coupling J. 

Site diluted disordered antiferromagnets have been studied previously in the frame- 
work of the NLSM, leading to very interesting results [Q. The aim of this work is to 
investigate the consequences of disorder in the quantum Heisenberg antiferromagnet 
by considering a continuous random distribution of spin stiffnesses ps in the NLSM. 
The effects of disorder are particularly interesting and, in fact, lead to unexpected 
results in the case of skyrmion correlation functions. Quantum skyrmion states | sk > 
are characterized by the property 

Q\ sk >= I sk > 

where Q is the topological charge operator. These excited states are created out of 
the ground state by an operator /i whose correlation functions and properties have 
been extensively studied in for the case of the NLSM. The skyrmion energy, in 
particular, can be inferred from the large distance behavior of the skyrmion correlation 
function < /i/i^ >. 

The pure NLSM at zero temperature is known to exist in two phases 0. These 
can be characterized by an order parameter, which is the ground state expectation 
value of the continuum limit of the sub lattice spin operator, namely < a >, along 
with a dual (disorder) parameter given by < p > 0. One of the phases of the pure 
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system is ordered, having < a >^ and < ^ >= 0. In this phase < fifi' > has an 
exponential large distance decay implying the existence of a nonzero creation energy 
for the skyrmions 0. The other phase is a paramagnetic quantum disordered one, 
presenting < a >= and < fi >y^ 0. In this phase, the fact that the skyrmion states 
|sk >= fi\0 > are not orthogonal to the vacuum mean that there are actually no 
genuine skyrmion excitation in the system. A third possibility, which is not realized 
in the pure NLSM at zero temperature, would be a phase in which < >= with 
the skyrmion correlation function presenting a power-law decay at large distances. In 
this case, the system would have zero energy skyrmions in its excitation spectrum 
and < a >= 0, the absence of order being closely related to the vanishing of the 
soliton energy. The possibility of both < fi >^ and < a >7^ 0, on the other hand, 
is forbidden by a duality relation existing between the spin and soliton operators, 
which has been rigorously demonstrated in one spatial dimension and, for physical 
reasons, should also be valid in higher dimensions. 

One can ask whether some different phases may occur when disorder is introduced 
in the system. Starting from the ordered phase of the pure NLSM, which corre- 
sponds to the Heisenberg antiferromagnet, we investigate this possibility by studying 
quenched averages in the presence of random couplings. When magnetic dilution is 
not exponentially suppressed, we conclude that the system exists in a new phase, 
where << /i >>= and << a (the quenched averages <<>> are defined in 

Section 3) but in which the skyrmion correlation function presents a power-law decay 
at large distances that implies the existence of zero energy skyrmions. This phase, 
which never occurs in the pure system, is a "soft" disordered Neel phase, possessing an 
order parameter << a >>^ 0, in spite of the fact that the skyrmion energy vanishes. 
The power-law decay of correlation functions in this phase shows some characteristics 
of criticality. Universality classes, which are determined by the behavior of the distri- 
bution of random couplings at p — > 0, also can be clearly identified in this new phase 
which bears some resemblance with a spin-glass phase. This is discussed in Sect. 6. 

When the distribution function of random couplings is such that magnetic dilu- 
tion is exponentially suppressed, on the other hand, the system is shown to exist in 
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a phase with « ji »— and << a »^ presenting, at the same time, an expo- 
nential decay of the skyrmion correlation function, analogously to what happens in 
the ordered Neel phase of the pure system. The skyrmion excitations always possess a 
nonzero energy in this phase. Only the subleading behavior of the skyrmion and spin 
correlation functions at large distances is modified by the introduction of disorder in 
this phase, which might be called a "hard" disordered Neel phase. 

The paper is organized as follows. In Section 2, we review some properties of the 
quantum NLSM relevant for the present work as well as the continuum limit of two- 
dimensional quantum antiferromagnets in the pure case. In Section 3, we consider 
disorder in the NLSM, manifested in a random distribution of couplings in the Neel 
phase. The probability distribution functions for these couplings are also introduced. 
In Section 4, we study the quenched averages of skyrmion correlation functions in the 
case where magnetic dilution is exponentially suppressed. In Section 5, we consider 
the same averages in situations in which magnetic dilution is not suppressed. We also 
show the occurrence of a new phase, presenting some characteristics of criticality, in 
which the system belongs to universahty classes determined by the behavior of the 
distribution function at ps 0. Discussion of the results, conclusions and future 
perspectives are presented in Section 6. 

2) The Quantum Nonlinear Sigma Model and the 
Pure Heisenberg Antiferromagnet 

2.1) The Quantum Nonlinear Sigma Model 

Let us start by reviewing the properties of the two-dimensional 0(3)-symmetric quan- 
tum NLSM. Subsequently, we shall recall how it is mapped in the 2D Heisenberg 
antiferromagnet. The NLSM is defined by the action [d^x — drdPz) 



where the field n is subject to the constraint = 1. po is a coupling parameter and c, 
a characteristic velocity. Henceforth, unless otherwise specified, we shall make fi—1 
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and c = 1. Writing the nonlinear sigma field as n = (cr, vf), the zero temperature 
partition function can be expressed as 



Z = J DaDnDXexpS^- J 



(fx 



-[{d,af + \d,n\ 



a 



Po 



(2.2) 



where we rescaled the fields and introduced the constraint trough the Lagrange multi- 
pher field A. We use the notation (9^ = (^, V). Integrating on vf, we get the effective 
partition function 

1 



(fx 



2(W + ^A 



cr -po 



DaDX exp 

+trln[-n+a]} (2.3) 
The constant saddle-point equations derived from the above expression are 



<a>' 



< X >< a >= 



Po 



(2.4) 



where ^ 



(27r)3 P + m2 

i < A >. At zero temperature, the system presents two phases 0: 
an ordered Neel phase, for which < A >= and < a >^ and a (paramagnetic) 
quantum disordered phase, in which < A >^ and < a >= 0. We explore below the 
physical properties of the basic excitations of the system and corresponding correlation 
functions in each of these two phases. 

An important feature of the NLSM is the existence of topological excitations, 
called skyrmions. Classically, they are solutions of the field equations carrying the 
topologically conserved charge 



Q = — d^xe'^e^^'n^diH^dX 



(2.5) 



At quantum level, the skyrmion states | sk > are eigenstates of the Q operator with 
eigenvalue equal to one and are created by an operator /i satisfying the commutation 
rule [Q, /i] = /i. The correlation functions of this operator have been studied in 
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detail in the ordered phase of the NLSM, taking into account full quantum effects 
0. Together with < a >, the ground state expectation value of the soliton creation 
operator ^ is a convenient tool for the characterization of the phases of the system, 
which we are going to exploit. 

2.2) The Quantum Disordered Phase 



We start with the quantum disordered phase, where m 7^ 0. Evaluating the integral 
in ( P?^ ) using the large k cutoff A and taking < a >= 0, it is easy to see that 

^ Po > (2.6) 



An 27r2 

Note that the large A behavior of ( p.6|) , as usual, my be compensated by the bare 
coupling pq, yielding a finite parameter m. Using the saddle-point solution z < A >= 

2 

it becomes clear that, up to a constant, the effective cr-field action is given by 



SeffW] = I dx- 



(d^ay + m^a^ (2.7) 



The spin correlation function, in this phase, therefore, is given by (notice that we are 
working with imaginary time r = it) 

dio r d^k (^k-x^iu T 



< a{x, t)(t(0, 0) >QD= — J 



(27r)2 ,^2+ 1^12 + ^2 



(2. 



47r[|f|2 + r2]i/2 

The exponential decay reveals the presence of a correlation length ^ = m~^. The 
large distance behavior < aa >qd — ^ confirms that < a >qd= in this phase. 
Conversely, no quantum soliton excitations are expected to present in this phase and 
therefore we must have < /i >^ 0, implying that the quantum skyrmion state | sk > 
is not orthogonal to the ground state. As we shall see this will be confirmed below. 
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2.3) The Ordered Phase 



We now turn to the ordered phase. In this case, we have m = 0. Evaluating the 
integral in ( p.4|) , again using the large k cutoff A, we get 

A 



< a > 



ORD- 



Po 



27r2 



> 



(2.9) 



Once more, the large A behavior in ( |2.9| ) can absorbed in a redefinition of the bare 
coupling pQ. Introducing the renormalized (finite) coupling ps, the spin stiffness, as 

A 



Po 



> 



(2.10) 



we see that < a >ord— Psi which is nonzero in the ordered phase. At the quantum 
critical point ps = 0, the system enters the disordered phase, where m ^ 0. The 
sublattice magnetization M is given by M =< a >. In the ordered phase, we have 
MoRD = y/Ps while in the quantum disordered phase, of course, Mqd = 0. 

Let us consider now the complete renormalization of the theory in the ordered 
phase. From ( |2.10| ), we can write 

A \ 



Po = Zps 



Z 



1 + 



Introducing the renormalized fields n^, Xji and action Sji through 

Z-'/^nR ; X = ZXn ; Sn = S + 6S 



n 



where 



5S = -i{Z-^ -1) J ZXr 
it is easy to see that the renormalized action is given by 



S 



R 



\^r\ 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



This is identical to the classical action, but with renormalized, physical quantities, 
replacing the bare ones. For this reason, this phase in known as "renormalized clas- 
sical" [^. We see, in particular, that the physical coupling constant in this phase is 
the spin stiffness Ps- 
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Replacing 5" for Sr in (|2.3|) , inserting the saddle-point value < Xr >= and 
shifting the u^-field around its vacuum expectation value < >= ^/p^, namely 
defining rj = — ^/pj we get 

Sefflv] = J d'x^id.vr (2.15) 

This is the well-known Goldstone boson action and the corresponding correlation 
functions are 

< r]{x,r)r]{0,0) >=< aR{x,T)aR{0,0) >oaf - < c^r >oaf~- 



47r[|f|2 + r2]i/2- 

(2.16) 

We now see that < aRaR >oaf — ^< o-r >oaf¥' 0; large distances, thus confirming 
the fact that < aR >oaf7^ in this phase. 

In the ordered antiferromagnetic phase, we have the occurrence of classical skyrmion 
excitations, possessing Q = 1. These are given by |^ 

ns{x) = ps (sin /(r)f, cos /(r)) (2.17) 

with 

/(r) = 2 arctan - 
r 

where / is an arbitrary scale and r is the radial distance in two-dimensional space. 
The energy of this classical skyrmion excitation in the ordered phase, described by 
the renormalized classical action (|2.14|) and measured with respect to the ordered 
antiferromagnetic background is E = Anps. This must compared with the full quan- 
tum result obtained from a quantized skyrmion field theory. The two-point quantum 
skyrmion correlation function has been evaluated in the ordered phase of the quantum 
NLSM § and the result is 

< p{x,r)p\0,0) >oAF= exp {-271 ps[\x\^ + T^f^} (2.18) 

From this we can infer that the actual energy of the full quantum skyrmions in the 
ordered antiferromagnetic phase is 27rps, that is, a half of the classical value. From 
( p.l8| ), we also see that < pp'^ >oaf — ^ 0, at large distances, implying that < p >= 
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in this phase. This means that the quantum skyrmion state | sk > is orthogonal to 
the vacuum and the quantum skyrmions are genuine excitations. For ^ 0, on the 
other hand, when we approach the quantum critical point leading to the disordered 
phase, we see from ( |2.18| ) that < /i >^ 0, confirming therefore our anticipation for 
the ground state expectation value of the skyrmion operator in the disordered phase. 

2.4) Connection with the Heisenberg Antiferromagnet 

Two-dimensional antiferromagnets on a square lattice can be described by the 0(3)- 
symmetric Heisenberg hamiltonian, given by 

^ = E -^ii^i ■ Sj (2.19) 

<ij> 

where the sum runs only over nearest neighbor sites and Jij > 0. The "pure" case is 
characterized by the fact that the coupling constants Jij are determined and fixed. In 
the homogeneous case, all the coupling constants are equal and we have Jij = J > 0. 
At zero temperature, this system is known to exist only in an ordered Neel phase. The 
quantum fluctuations are not capable of destroying the long-range antiferromagnetic 
order for any value of the coupling constant [0]. 

It has been shown that in the continuum limit, the above quantum hamiltonian, 
in the homogeneous case, is mapped into the ordered phase of the quantum NLSM 

0, the nonlinear sigma field n(x, t) being the continuum limit of the sublattice spin 
operator. The spin stiffness ps which, as we saw, controls the whole physical properties 
of the system in the ordered phase, is related to the Heisenberg antiferromagnetic 
coupling J as 0, ^ 

Ps = JS^Zp^ (2.20) 

where S is the spin quantum number and Zp^ is a constant accounting for quantum 
corrections to the classical continuum limit. For 5* = |, we have [0, |[ 

~0.18J (2.21) 

In the non-homogeneous case, where the coupling constants Jij are different for 
each link, we can derive the continuum limit by following the same procedure as in 
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provided the configuration of coupling constants Jij in ( p.l9| ) is slowly varying (this 
is going to be made precise in what follows). In this case, we obtain in the continuum 
limit, a Nonlinear Sigma Model with the spin stiffness ps replaced by a slowly varying 
configuration p{r) that is related to in the same way that ps is related to J, namely 

S = Jd'x [l\d,n\' + i\ [|n|2 - p,(r)\ } (2.22) 

This is equivalent to modifying the constraint from (5[|np — 1] to (^[|np — /(r)], with 
Poif) = pof{r). Integrating over tt, we get 

S = J d^x!^^{d^af + i\[a^ - po{f^]^ +tr\n[-n + tX\ (2.23) 

Now translation invarianace is lost and the saddle-point equations become 

-V^ <a> +2m^{r) < a >= 

<->^=P«W-/(S5JI^-^>W (2.24) 

We see that now < a >=< a > (r) and consequently m{r) ^ 0, in spite of the fact 
that < o 0. The spin correlation function becomes damped now, with a damping 
factor m{r). This is in agreement with previous investigations of spin- waves in similar 
situations ^ 

The soliton correlation function can be evaluated in a saddle-point approximation 
in the theory described by (|2.22|) [^]. This has a simple expression in terms of 
< a >, which will be very convenient for the obtainment of quenched averages in the 
disordered version of the model, namely, 

< /i(f , t)p\0, 0) >= exp [-2n < a >^ [|f ^ + r^^^} = exp [-27rp{f)[\x\^ + r^^^} 

(2.25) 

In the rest of this work, we consider the situation in which intrinsic disorder is 
introduced in the system and investigate its effects on the soliton correlation function. 

3) The Continuum Limit of Disordered Antiferro- 
magnets 
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3.1) The Disordered System 



Let us describe the presence of disorder in the two-dimensional Heisenberg quantum 
antiferromagnet, by considering a random distribution of couphngs Jjj in ( ^.191 ), in 



analogous way as in the Edwards- Anderson model [|Tl|]. Here, however, we will keep 
only antiferromagnetic couplings Jij > 0. This will allow us to easily obtain a con- 
tinuum field theory version for the disordered model, in the same way as in the pure 
case. The disorder is introduced in the continuum version by taking a random dis- 



tribution P[p{f)] for the slowly varying spin stiffness p{r) appearing in ( p.22| ). We 



require that P[p{r)] = for p{r) < 0. This will ensure that in spite of the presence of 
disorder, we are always in the ordered phase of the NLSM, for which the mapping to 
the Heisenberg antiferromagnet exists. We also impose the condition that the vari- 
ance of this distribution is always much smaller than ps, in order to ensure that the 
p(r)-configurations are slowly varying. 

We are only going to consider the quenched case and take quantum averages at 
zero temperature using a fixed configuration for p{r). Subsequently we shall evaluate 
the average over the p(r)-configurations using the P[p{r)] distribution function. The 
relevant average for an operator A in the quenched random system will be therefore 

« A jUr^ dp{r)P[p{r)] <A>, (/9(r-)) (3.1) 

where A = (H^- 1) is a normalization factor corresponding, in the lattice, to a product 
over all links {ij) and < ^ >g is the zero temperature quantum average. 

3.2) Distribution Functions 



Let us introduce now the distribution functions we are going to use, in order to 
describe the disorder. We shall consider basically two functions containing a gaussian 
distribution centered around the pure spin stiffness ps- The first one is (henceforth 
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we omit the argument r in p) 

(3.2) 



p<0 
where u > 0. In this expression, we assume that both A << ps and 6 << ps, thereby 
guaranteeing that the random p-configurations are slowly varying. We also assume 
that 

where L is the maximum dimension of the system and c, the spin-wave velocity, its 
characteristic velocity. Two regimes of disorder described by (|3.2| ) can be distinguished 
and will produce a completely different behavior of the correlation functions, as we 
shall see. A first one is obtained by choosing a symmetric gaussian, with 6 = A. A 
second one with the choice S « A and 

In the second case, there is a severe exponential suppression of values of the spin 
stiffness around ps = 0, that is, magnetic dilution is exponentially suppressed in a 
very strong way. In the first case, dilution is not so much suppressed and, as we shall 
see, the system has the same qualitative behavior as a diluted one. 

Experimental values for the parameters of the above distribution, in the case of 
high-temperature superconducting cuprates in the ordered antiferromagnetic phase, 
which are typical examples of two-dimensional Heisenberg antiferromagnets are 
Ps — lO^^eV, he ~ leVA. For a sample of dimension L ^ 1mm, we choose A ~ 
10~^eV, which satisfies ( p.3|) . For 6, we have 5 = A in the case with dilution. In 
the case without dilution the choice 6 ~ lO^^eV will satisfy (|3.4|) . In both cases, the 
condition A,S << ps is satisfied. 

In the distribution function ( p.2| ), the normalization factor is given by 

Ni = (^0 [A'^ + 6"] (3.5) 

12 



r(^) 
P = Ps+ } I { A 



The average spin stiffness is 



in tlie case wliere dilution is not suppressed and 



(3.6) 



P = P. 



+ V2- 



r(| 



A 



(3.7) 



in the presence of exponential suppression of dilution {6 << A). 
The second distribution function we are going to use is 

P > 



N2 



2^ 



(3.8) 



p < 

where u > 0. We assume a « ps, again to ensure that the p-configurations are 
slowly varying. This condition can be experimentally satisfied, in the case of the 
high-temperature cuprates, with a choice of cr ~ lO^^eV. Observe also that using 
the same experimental values of the previous paragraph, this value of a satisfies the 
condition 

which is similar to ( ^.31 ). 

The normalization factor in ( |3.9| ) is now given by 

2 

N2 = a'^r(j.)D_,(-^)e-^ (3.10) 
a 

where D^u^x) is a parabolic cylinder function. For the P2[p] distribution function, 
the average spin stiffness is 

p = p, + (z/-l)— (3.11) 

Ps 



4) Disorder With Exponentially Suppressed Dilu- 
tion 



In this section, we consider the situation in which magnetic dilution is exponentially 
suppressed in the disordered system. As explained above, this correponds to the choice 
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of Pi[p] as the distribution function , with the parameters A and S satisfying (|3l3| ) 
and ( p.4|) , respectively. We are going to evaluate the quenched skyrmion correlation 
functions, starting from the ordered antiferromagnetic phase of the pure system. In 
this case, the zero temperature pure quantum averages < /i/i^ >oaf, given by (|2.18D, 



depend exponentially on the spin stiffness. For the disordered NLSM, introduced in 
subsection 2.2, in the regime where the spin stiffness p{r) is slowly varying, the soliton 



correlation fubction is given by ( |2.25| ). Hence, when evaluating the quenched averages 
we shall have a zero temperature quantum average whose p(r)-dependence is of 
the form exp{— ap(r)} where a = 27iX, with X = + r^]-*^/^. The relevent integral 
for the evaluation of ( |3.1| ) is, therefore 

A = Ai + A2 (4.1) 

where 



A = — / dpips-py-'e-'^^e-^^ = -— dxiaps-xy-'e-^e-^^ 
JMi Jo JMia" Jo 



(4.2) 



and 



^ = — / dpip - p^y-'e-'^'^e-^^ = -— / dxx^-^e-^e-S^ (4.3) 
i\i Jps JMia" Jo 

We shall be interested in the behavior of quenched averages of correlation functions, 

given by ( ^.1[ ), at large distances {X —>■ oo). In this case we have a —>■ oo and we can, 

therefore, use conditions (|3.4|) and (p.3|) , respectively, in ( [4.3|) and ( |4.2|) , to obtain 

and 

Nia'' Jo Ni a" ^ ^ 

(we could also have obtained the result in (|4.5|) by exactly evaluating the last integral 

in (^^ ) and subsequently taking the limit a — » (see ( ^.7|) and ( |5.8| ))). We conclude 

from and (U) that 

+ — 



p-(^Ps 



(4g) 
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where we used (|3.4|) in the second hne, in order to get the dominant behavior for large 
a. 



tance behavior of the skyrmion quenched averages in the absence of magnetic dilution, 
namely, 



Observe that the presence of disorder, in the case where dilution is esponentially 
suppressed, does not modify the dominant exponential large distance behavior. Only 
the subdominant power law decay is modified by the introduction of an additional 
exponent v. 

5) Disorder Without Exponential Suppression of 
Dilution 

In this section, we are going to consider the situation in which magnetic dilution is 
not exponentially suppressed by the random distribution of couplings. As we saw in 
Sec. 2, this can happen either when we use the distribution Pi{p\ in the symmetric 
case when 5 = A or when we use P2[p]- In what follows, we study the two cases 
separately. 

5.1) Distribution Function Pi[p\ 

In this case, the relevent integral for the evaluation of ( |3.1| ) is 



From ( p.l8| ) and (O) we can immediately obtain the expression for the long dis- 




(4.7) 



B = Bi + A2 



(5.1) 



where A2 is given by ( |4.3| ) and 



Bi = — dp{ps - py-'e-'^Pe-^^ = —— / dxx'^-^e^e-^ 




(5.2) 



In the large distance regime, when a 00, we can use (|3.3|) to obtain [0 




-ap: 



.V 



-ap: 



(5.3) 
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where iFii^u; 1 + z/; aps) is a confluent hypergeometric function. Using the large dis- 



tance asymptotic behavior of this function, we get [p!2 



Po 



Nia 



Combining (5^) with ( [4.6|) , we see that 



B 



Ps 



-apa 



a—*oo Hs 



y-1 



Nia 



(5.4) 



(5.5) 



a a" 

From this result, we can immediately infer the large distance behavior of the skyrmion 
quenched correlation function. This is given by 



t X^oo Ps 

« /i/i' »OAF - 



2m \x) ■ <5-«) 

Now the large distance behavior of the correlation function is drastically changed. 
The previously dominating exponential decay is completely washed out and we have, 
instead, a power-law decay. As we shall argue in the next subsection, the exponent of 
the power-law is universally determined by the behavior of the distribution function 
at p — > 0. 

5.2) Distribution Function P2[p] 



Let us consider now the situation in which the disorder is described by the distribution 

is now 



function P2[p]. The relevant integral for the evaluation of (|0 



12 



1\2 Jo 



X exp 



a 



ps 



exp 



D^i, aa 



(5.7) 



where D^^{x) is the parabolic cylinder function. The large distance behavior of (|5.7| ) 
can be obtained by considering the property 



(5. 



D_,,(—x] — ^ „ , , e 4 X 
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Using this and (|3.9| ), we can immediately obtain the large distance behavior of the 
quenched skyrmion correlation functions for the distribution function P2[p\- This is 
given by 



We observe here that, as in the case of the distribution used in the previous subsec- 
tion, the introduction of disorder completely modifies the large distance behavior of 
the correlation functions, eliminating the exponential decay. This fact can be gen- 
erally understood by observing that the large distance (large a) behavior of (|3.1|) is 
determined by the behavior of the distribuiton P[p] at p ^ 0. This happens because, 
for the distribution functions P[p\ (with support in the region p > 0) and quantum 
averages < A >q (p) (exponentially depending on p) considered in this work, the 
quenched average ( |3.1[ ) is proportional to the a Laplace transform of P[p\- As a con- 
sequence, the asymptotic large distance behavior of (|3.1| ) is universally determined by 
the behavior of P[p] for p — 0. We see, for instance, that the values of the parameter 
V in P2 [p] determine universality classes to which the disordered system belongs. This 
can be confirmed by the behavior of the correlator ( ^.6|) . Also, by comparing (|5.6| ) 
with ( |5.9| ) we conclude that distribution Fi[p], in the symmetric case when 5 = A, 
is in the p = 1 universality class of P2[p\- This can be also verified by observing 
that both distributions have the same type of behavior at p ^ 0. Conversely, in 
the asymmetric case when 5 << A, the distribution function Pi[p\ is exponentially 
suppressed at p ^ 0, a type of behavior which is never presented by P2[p], leading to 
an exponential decay of the quenched correlators. 

6) Discussion and Conclusions 

Our continuum analysis of two-dimensional quantum antiferromagnets in the pres- 
ence of a random distribution of couplings at zero temperature, has shown that the 
quenched averages of soliton (skyrmion) correlation functions are modified with re- 
spect to the pure case. The modification is particularly drastic whenever magnetic 
dilution is not exponentially suppressed in the disordered system. This effect has quite 




(5.9) 
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interesting consequences in the physical properties of skyrmions. It is a well known 
fact that an exponential decay of the soliton correlation function at large distances 
would indicate that the energy of the soliton excitations is nonzero and proportional 
to the coefficient of the exponent. In the pure NLSM in the ordered phase, this is 
given by Es — 27rps = 27r < cr A nonzero soliton energy, therefore, is associated 
to an ordered ground state with < a >^ 0. Physically this can be understood as 
a consequence of the fact that in an ordered ground state there is an energy cost to 
make the spin flips necessary for the introduction of a soliton state. The exponential 
decay of < /i/i"'' >, further implies through < fifj.'^ >— >• | < > p that < /v, >= 
which means that the soliton states are orthogonal to the vacuum, that is to say, true 
excitations. A power-law decay, on the other hand, while still leading to < fj, >— 
and therefore meaning that quantum solitons are genuine excitations would imply 
that the energy necessary for the creation of these solitons is equal to zero. In a 
generic pure system at zero temperature, this would correspond to a quantum disor- 
dered ground state (< a >— 0), because when the ground state is not an ordered 
one, there is no energy cost for introducing the spin flips necessary to create a soliton 
state. These zero energy quantum skyrmions, in spite of bearing few relation with 
their classical ancestors occurring in an ordered phase, would exist as true physical 
excitations in such a phase. It should be stressed, however, that the kind of quantum 
disorder occuring when we have a power-law decay of the soliton correlation function, 
differs from the one found in a paramagnetic phase, such as the quantum disordered 
phase of the NLSM, in which we have < ////^ >^ C 7^ 0, that is < fi >^ and 
< a >— 0. Consequently, we conclude that a power-law decay of the soliton corre- 
lation function would imply some different type of disorder than the one found in a 
paramagnetic phase. In summary we have the following possibilities for the phases of 
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a pure quantum antiferromagnet: 
< ^fJ >^^—^ C 7^ ; < ^ >7^ 0, < (J >= — Paramagnetic Quantum Disordered 

< jjfi^ >^^—^ e~^^^ ; < /i >= 0, < a — Antiferromagnetic (Neel) 

< ^jJ >^^~^ Yu ] < 1^ >= 0, < (J >= — Non Paramagnetic Quantum Disordered 

(6.1) 

The first two ones are realized in the pure NLSM at zero temperature and only the 
second one in the pure two-dimensional Heisenberg quantum antiferromagnet with 
nearest neighbors interaction, also at T = 0. 

We then consider the presence of disorder. Starting from the ordered antiferromag- 
netic phase of the pure NLSM, which corresponds to the Heisenberg antiferromagnet, 
we have studied in this work the effects of disorder introduced through a continuum 
random distribution of couplings. In all cases considered here, the quenched magne- 
tization is nonzero, namely, Mq =<< a >>~ ^/p, where p is given by ( |3.(j| ) ( PT^ ) 
and ( |3.11|) , respectively. The studied systems are always in an antiferromagnetic "dis- 
ordered" Neel phase having a nonzero order parameter in spite of the fact that the 
couplings are random. The types of disorder considered here, ab initio cannot destroy 
the antiferromagnetic order of the pure system since they only allow the presence of 
positive or null couplings. There are, however, two possible types of such phases, 
which we call "soft" or "hard", according to whether magnetic dilution is exponen- 
tially suppressed or not. 

When dilution is exponentialy suppressed, only the subleading term of the cor- 
relation functions at large distances is modified by disorder. The exponential decay 
of the pure system is preserved and the skyrmion energy in the disordered system 
is the same as in the pure case. The system exists in a phase corresponding to the 
second possibility in ( |6.1| ), which may be called a "hard" disordered Neel phase. Con- 
versely, when dilution is not exponentially suppressed, the large distance behavior is 
drastically changed from an exponential to a power-law decay. The skyrmion energy 
consequently becomes zero in the disordered system, despite the fact that the mag- 
netization is non- vanishing and << a »^ 0. This is a completely novel behavior 
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for quantum skyrmions. Physically, we can understand it as follows: the disorder 
introduced by the random distribution of couplings is sufficient to reduce to zero the 
energy necessary for the creation of a skyrmion configuration, even though it can- 
not destroy the quenched magnetization. The soliton correlation functions have a 
universal behavior, characterized by universality classes which are determined by the 
behavior of the distribution function at ^ 0. 

It would be interesting to find a physical realization for this phase in which 
skyrmion correlators have a power-law decay at large distances and which is associ- 
ated to a diluted disordered antiferro magnet. One might be tempted to associate it to 
a spin-glass phase since spin-glasses can be obtained by dilution of antiferromagnetic 
systems. A well known example in three dimensions is the substance ZnCr2-xGax04, 
obtained by diluting the pure antiferromagnet ZnCr20i with the nonmagnetic Ga 



atoms |TJ, 0. The possibility of the "soft" phase, obtained in the present work by 
dilution of a pure antiferromagnet, being a spin-glass, however, is ruled out by the fact 
that the quenched magnetization is nonvanishing. A key ingredient for this would be 
the presence of frustration which is absent in the system studied here but, conversely, 
present in the case of ZnCr2^xGaxOi. The diluted phase studied here is actually a 



new one which could be added to the list (|6.1| ) and is characterized by 



< nii^ ^ < ^ >= 0, < a >7^ - Soft Disordered Neel (6.2) 

In this phase, the skyrmion energy vanishes, not as a consequence of quantum fiuctu- 
ations, as in the third phase listed in ( |6.1| ), but rather, because of the type of disorder 
introduced when magnetic dilution is not exponentially suppressed. This, in spite of 
not being capable of destroying the order parameter, forces the skyrmion energy to 
vanish. The physical properties of the zero energy quantum solitons occuring in these 
kind of phases may play an important role in planar antiferromagnetic systems such 
as high temperature superconducting cuprates an shall be exploited elsewhere. For 
this purpose, one should investigate the effects of disorder in the continuum models 
for doped antiferromagnetic planar systems used to describe materials such as LSCO 
and YBCO |T|. 
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Let us stress, finally, that that the kind of magnetic dilution considered in this 
work, even in the cases where it is not exponentially suppressed, is always softer than 
a delta function-type of dilution containing a piece P[p] = xd{p), which would be suit- 
able for describing the type of dilution occuring in compounds such as La2Cui_j;ZnxOi 
@]. As future extensions of this work, we intend to consider more general disorder 
distributions including these and ferromagnetic couplings as well. In the latter case, 
we can expect to describe spin-glass phases as well. It is likely that in a spin-glass 
phase the behavior of soliton correlation functions shall be similar to the one found 
in the soft Neel phase of our model. We are presently investigating this point. 
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